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Abstract 

A  relationship  between  a  flooding  and  a  cell  voltage  drop  for  polymer  electrolyte  fuel  cell  was  investigated  experimentally  and  numerically.  A 
visualization  cell,  which  has  single  straight  gas  flow  channel  (GFC)  and  observation  window,  was  fabricated  to  visualize  the  flooding  in  GFC.  We 
ran  the  cell  with  changing  operation  condition,  and  measured  the  time  evolution  of  cell  voltage  and  took  the  images  of  cathode  GFC.  Considering 
the  operation  condition,  we  executed  a  developed  numerical  simulation,  which  is  based  on  multiphase  mixture  model  with  a  formulation  on  water 
transport  through  the  surface  of  polymer  electrolyte  membrane  and  the  interface  of  gas  diffusion  layer/GFC.  As  a  result  in  experiment,  we  found 
that  the  cell  voltage  decreased  with  time  and  this  decrease  was  accelerated  by  larger  current  and  smaller  air  flow  rate.  Our  simulation  succeeded  to 
demonstrate  this  trend  of  cell  voltage.  In  experiment,  we  also  found  that  the  water  flushing  in  GFC  caused  an  immediate  voltage  change,  resulting 
in  voltage  recovery  or  electricity  generation  stop.  Although  our  simulation  could  not  replicate  this  immediate  voltage  change,  the  supersaturated 
area  obtained  by  our  simulation  well  corresponded  to  fogging  area  appeared  on  the  window  surface  in  the  GFC. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Water  management,  to  control  the  water  transport  to  avoid 
dehydration  and  flooding,  is  critical  issue  to  improve  the  perfor¬ 
mance  of  polymer  electrolyte  fuel  cell  (PEFC).  Operating  PEFC 
under  high  current  density  and/or  low  Stoichiometric  condition 
often  invites  the  flooding,  which  forms  locally  and  spreads  in 
the  cell.  Once  the  flooding  starts,  the  reactants  gas  supply  to  cat¬ 
alyst  layer  is  blocked,  and  then  the  cell  performance  deteriorates 
accompanying  with  increase  of  concentration  overpotential.  To 
improve  the  performance  of  PEFC,  understanding  the  mech¬ 
anism  of  the  water  transport  in  cell  is  required  as  first  step, 
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and  then  the  proper  design  and  the  optimization  of  operation 
condition  should  be  established. 

The  flooding  appears  in  every  segment  in  cell,  such  as  gas 
flow  channel  (GFC),  gas  diffusion  layer  (GDL)  and  catalyst  layer 
(CL).  The  flooding  in  GFC  was  observed  with  visualization  cell, 
whose  separator  partially  made  with  transparent  material  such 
as  glass  [1,2].  In  addition  to  the  visualization  in  GFC,  recent 
neutron  imaging  technique  [3]  enables  us  to  visualize  whole 
region  in  cell.  Though  it  is  generally  difficult  to  visualize  the 
flooding  in  GDL,  which  is  characterized  as  micro-porous  and 
opaque  material,  neutron  beam  property  and  advanced  spec¬ 
troscopy  technology  enable  us  to  visualize  every  segment  in 
cell. 

The  numerical  simulation  of  flooding  in  PEFC  has  been  also 
progressed.  Based  on  the  multiphase  mixture  model  [4],  a  lot  of 
two-phase  flow  simulations  of  PEFC  have  been  reported:  a  sim¬ 
ulation  which  turns  its  attention  especially  to  cathode  channel 
and  GDL  [5]:  a  simulation  explicitly  considering  catalytic  layer 
[6] :  a  simulation  focused  on  the  GDL  under  rib  [7] :  a  simulation 
for  transient  response  of  unit-cell  [8]. 
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Nomenclature 

A  area  (m2) 

c  concentration  (mol  m- 3 ) 

C  capacitance  (F) 

D  diffusion  coefficient  (m2  s  - 1 ) 

F  Faraday’s  constant  (C  mol-) 
i  current  density  (A  cm-2) 

I  current  (A) 

k  permeability  (m2) 

K  transport  coefficient  (ms-1) 

ky  phase  change  rate  constant  (1  s-1) 
mg i  phase  change  rate  (kg  (m3  s)-1) 

n&  electro-osmotic  coefficient 

P  pressure  (Pa) 

Q  electric  charge  (C) 

R  gas  constant  or  resistance  (J  (mol  K)-1  or  Q) 

s  liquid  water  saturation 

T  temperature  (K) 

u  velocity  vector  (ms-1) 

V  volume  or  cell  voltage  (m3  or  V) 

Y  mass  fraction 

Greek  symbols 
8  thickness  (pan) 

s  porosity 

rj  overpotential  (V) 

A  relative  mobility 

/x  viscosity  (Pa  s) 

v  kinematic  viscosity  (m2  s-) 

0C  contact  angle  (°) 

p  density  (kg-m-3) 

o  surface  tension  or  proton  conductivity  (N  m-1  or 

Sm-1) 

t  tortuosity 

§  utilization  ratio 

Subscripts 

ano  anode 

act  active 

BV  Butler-Volmer 

c  capillary 

cat  cathode 

eff  effective 

eq  equilibrium  state 

g  gas  phase 

gdl  gas  diffusion  layer 

gfc  gas  flow  channel 

1  liquid  phase 

m  membrane 

OC  open  circuit 

ohm  Ohmic 

rep  representative 

sat  saturation 


Superscripts 

h2 

hydrogen 

h2o 

water 

I 

gas  species 

liq 

liquid 

o2 

oxygen 

P 

pore 

vap 

vapor 

w 

water 

Thus  the  both  of  visualization  experiment  and  numerical  sim¬ 
ulation  for  the  flooding  have  been  progressed.  However,  precise 
cross  comparison  between  them  does  not  seem  to  have  been 
done.  Almost  numerical  studies  have  been  evaluated  only  by  the 
comparison  with  the  measured  current-voltage  characteristics, 
which  is  categorized  as  whole  and  steady  state  characteristics. 
Liquid  water  saturation  and  its  time  evolution  has  not  been  com¬ 
pared  between  simulation  and  visualization  study. 

On  the  basis  of  this  background,  a  simple  visualization  cell, 
which  has  single  straight  GFC,  was  fabricated  to  precisely  com¬ 
pare  numerical  result  with  the  visualized  flooding  images  in 
GFC.  We  ran  the  cell  with  changing  operation  condition,  and 
measured  the  time  evolution  of  cell  voltage  in  addition  to  suc¬ 
cessively  recording  the  visualized  image  in  GFC.  Considering 
the  geometry  of  the  cell  and  the  operation  condition,  we  executed 
a  developed  numerical  simulation,  which  is  based  on  Multiphase 
mixture  model  with  the  punctilious  formulation  on  water  trans¬ 
port  through  the  polymer  electrolyte  membrane  (PEM)/GDL 
and  the  GDL/GFC  interface.  And  then  we  evaluated  the  numer¬ 
ical  simulation  result  by  the  comparison  with  the  visualization 
experiment  result.  In  the  following,  the  visualization  experi¬ 
ment  method  and  the  formulation  used  in  the  simulation  are 
explained.  Subsequently  experimental  result  and  numerical  sim¬ 
ulation  result  are  shown. 

2.  Experimental  method 

2.7.  Cell  specification 

To  investigate  the  liquid  water  behavior  in  GFC,  a  simple 
visualization  cell,  which  has  single  straight  GFC,  was  fabricated 
as  shown  in  Fig.  1.  Comparing  typical  cell  such  as  serpentine 
and  parallel  flow  channel  cell,  our  simple  cell  has  advantages:  we 
can  directly  analyze  the  relationship  between  flooding  and  cell 
voltage  without  additional  consideration  such  as  water  transport 
under  rib  appeared  in  typical  cell:  the  cell  geometry  can  be  easily 
incorporated  into  numerical  simulation. 

Fig.  1  shows  top  and  side  view  of  the  visualization  cell.  Sepa¬ 
rators  and  back-plates  were  made  of  SUS3 16,  and  the  separators 
were  gilded  to  reduce  electric  contact  resistance  and  to  protect 
separators  from  corrosion.  A  rectangular  GFC  was  formed,  and 
its  width,  depth  and  length  are  1.6,  1.0  and  30  mm,  respectively. 
The  geometry  and  the  size  of  anode  channel  are  as  same  as  that  of 
cathode  channel.  We  used  Gore  PRIME  A®  membrane-electrode 
assembly  (Series  5510,  electrolyte  membrane  thickness:  30  [xm, 
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Back  plate 


©Gas  flow  channel  (GFC) 

Separator  (Shoulder)^ 


1 


V 


Inlet  (6) 


Gas  diffusion  layer  (GDL) 


*(5)  Outlet 


MEA  Gasket  Separator 


Fig.  1.  Top  and  side  view  of  visualization  cell. 


Table  1 

Dimension  of  visualization  cell 


Quantity 

Value 

Width  of  GFC  (mm) 

1.6 

Depth  of  GFC  (mm) 

1 

Length  of  GFC  (mm) 

30 

Width  of  GDL  (mm) 

3.2 

Thickness  of  GDL  (pan) 

300 

Length  of  GDL  (mm) 

30 

Active  area  (cm2) 

0.96 

Thickness  of  PEM  (pm) 

30 

catalyst  loading  of  anode  and  cathode:  0.3  mg  cm-2),  which  was 
sandwiched  between  two  sheets  of  GDL  (CARBEL-CL®,  car¬ 
bon  cloth  type).  GDL  was  cut  as  3.2  mm  x  30mm  so  that  it 
covers  GFC  area.  We  assume  that  this  GDL  area  of  0.96  cm2 
determines  the  effective  active  area.  Other  dimensions  of  the 
cell  are  listed  in  Table  1 . 


2.2.  Test  bench  for  cell 

Fig.  2  is  the  cell  test  bench,  which  consists  of  the  control  sys¬ 
tem  of  gases  flow,  temperature  and  electricity.  The  reactant  gas 
for  cathode  was  air,  and  that  for  anode  was  hydrogen.  Mass  flow 
controller  (CMQ-V,  YAMATAKE)  regulated  gases  flow  rate. 
Both  anode  and  cathode  gas  were  humidified  with  the  bubbler, 
whose  temperature  was  controlled  by  water  bath.  The  separa¬ 
tor  and  the  ambient  of  cell  were  kept  at  constant  temperature. 
The  visualized  image  in  GFC  was  recorded  through  transparent 
window  of  cell  with  microscope  (VHX-200,  KEYENCE).  AC 
impedance  (1280B,  Solatron)  was  also  recorded  intermittently 
to  grasp  Ohmic  and  activation  overpotential. 


gas  port 


mass  flow 
controller 


humidifier 

Fig.  2.  Test  bench  for  cell  operation. 


hot  air  electric  load 


Table  2 


Operation  conditions 


Number 

Current  (A) 

Air  flow  (seem) 

H2  flow  (seem) 

Utilization  ratio 

1-1 

0.15 

58 

22 

0.05 

1-2 

0.15 

116 

46 

0.025 

1-3 

0.15 

232 

92 

0.0125 

2-1 

0.3 

58 

22 

0.1 

2-2 

0.3 

116 

46 

0.05 

2-3 

0.3 

232 

92 

0.025 

3-1 

0.6 

58 

22 

0.2 

3-2 

0.6 

116 

46 

0.1 

3-3 

0.6 

232 

92 

0.05 

Cell  temperature  (°C):  30,  inlet  relative  humidity  (%):  75. 


2.3.  Operation  conditions 

A  main  purpose  of  our  experiment  is  to  observe  flooding 
in  GFC  under  a  given  operation  condition,  and  then  try  to 
understand  the  impact  of  the  flooding  on  cell  voltage.  The  cell 
temperature  was  controlled  to  be  30  °C,  which  is  considered  to 
be  rather  low.  This  low  cell  temperature,  however,  allows  us 
to  observe  flooding  phenomena  efficiently,  because  lower  tem¬ 
perature  promotes  to  condense  water  vapor  produced  by  cathode 
reaction.  The  relative  humidity  of  reactant  gases  at  GFC  inlet  was 
controlled  to  be  75%.  Other  operation  conditions  such  as  load 
current  density  and  flow  rate  are  listed  in  Table  2.  Our  simple 
cell  was  somewhat  unstable.  So  we  chose  rather  high  utiliza¬ 
tion  ratio  in  this  experiment.  The  visualized  image  in  GFC  with 
microscope  was  recorded  at  every  15  s. 

Comparing  with  practical  operation  condition,  the  cell  tem¬ 
perature  and  the  utilization  ratio  in  this  study  are  low  and  high, 
respectively.  However,  it  is  thought  that  our  result  is  still  helpful 
for  the  practical  case.  The  water  behavior,  which  is  shown  in  the 
later  section  and  denoted  as  flooding  and  flushing,  also  appears 
under  practical  condition.  Water  transport  in  our  case  is  thought 
to  be  essentially  same  as  practical  case.  All  we  have  to  pay 
attention  is  the  nature  that  lower  temperature  and  the  smaller  uti¬ 
lization  work  to  accelerate  water  production  and  water  drainage, 
respectively. 

3.  Mathematical  model 

The  simulation  domain  is  shown  in  Fig.  3.  A  unit-cell  with 
single  straight  GFC  is  numerically  simulated  as  two-dimensional 
field.  The  simulation  domain  includes  GFC,  GDL  and  PEM. 
The  multiphase  mixture  model  [4]  was  used  for  two-phase  flow 
formulation  in  GDL,  and  Stefan-Maxwell  equation  [9]  is  used 
for  multi-component  diffusion.  The  equivalent  electric  circuit 
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Fig.  3.  Calculation  domain.  Mass,  momentum  and  electric  transport  is  treated 
to  be  two-dimensional. 


Table  3 


Source  terms 

in  governing  equations  in  GFC  and  GDL 

Su 

Si 

Channel 

0 

—s*mg i  (i  =  vapor) 

GDL 

pe*u 

0(i  =  H2,O2) 

Here,  p  is  the  density  of  two-phase  mixture  defined  as, 

p  =  (\-s)pg  +  sp\  (3) 

Momentum  conservation  equation  for  the  mixture  in  GDL 
and  GFC  is, 


[10]  shown  in  Fig.  4  was  employed  to  calculate  the  cell  voltage, 
overpotentials  and  current  distribution.  Our  simulation  used  the 
following  assumptions:  (i)  CL  is  assumed  to  be  infinitely  thin 
and  the  cell  reaction  occurs  at  the  interface  between  PEM  and 
GDL,  (ii)  cell  temperature  is  constant  and  uniform,  (iii)  transport 
property  of  GDL  is  isotropic. 

Simulation  with  the  assumption  about  constant  and  uniform 
temperature  might  lead  to  the  overestimation  of  water  saturation. 
Especially  for  higher  current  density  condition,  this  overestima¬ 
tion  of  water  saturation  is  thought  to  be  promoted  considering 
that  the  high  current  increases  temperature  and  saturation  vapor 
pressure.  However  qualitative  impact  of  load  current  and  flow 
rate  on  water  behavior  can  be  simulated  at  least,  even  if  the 
constant  and  uniform  temperature  are  assumed.  This  point  is 
reinforced  by  the  correspondence  between  simulation  and  exper¬ 
iment  as  shown  in  the  later  section. 

3.1.  Governing  equations 

The  governing  equations  used  here  followed  our  previous 
study  [8]  and  the  reference  of  [4,9,10].  In  the  equations,  liquid 
water  saturation  of  s  was  used  as  key  variable,  which  is  defined 
as  the  ratio  of  liquid  water  volume  to  pore  volume: 

V\ 

s=—  (1) 

Vp 

Continuity  equation  for  gas  and  liquid  mixture  in  GDL  and 
GFC  is, 

^(e*p)  +  V-(pu)  =  0  (2) 

ot 


d 

dt 


+  v- 


( i  v 

VP  +  V  • 

(4) 


The  last  term  of  RHS,  which  can  be  treated  as  source  term, 
means  Darcy  flow.  In  GDL  and  GFC  region,  mass  conservation 
equation  of  gas  species  is, 

9 

-  [e*(l  -  s)pj]  =  -V  •  [(1  -  s)pgyl 'ug]  -  Vjg  +  St  (5) 
and  that  of  liquid  phase  is, 

^[e*pi£]  +  V  •  [puns]  =  s*mg\  (6) 

Mass  conservation  for  water  in  PEM  is, 

r)nw  7 

=  V  ■  (DmVpZ)  +  nd-MH2°  (7) 

ot  r 

The  effective  porosity  of  GDL  denoted  as  s*(=£/z)  used  in  above 
equations  becomes  unity  in  case  of  GFC.  The  source  terms  of  S 
shown  in  Eqs.  (4)  and  (5)  are  listed  in  Table  3. 


3.2.  Constitutive  relations 

3.2.1.  Diffusive  flux  of  multi- component 

Diffusive  flux  of  multi-component  Jg  in  Eq.  (5)  is  expressed 
as  [9], 


X 


=  -Ce 


d\\  d\2 
dix  d22 


Vx1 

Vx2 


d\\ 


x'D22 +x2D'3+x3D’2, 
Dll(VD%  +  (l-x')D'h 


dn  = 


di\  = 


r1  D23ID13  —  D12! 


(8) 


fi 


r2  D13tD23  —  DiZl 


12\ 


di2  = 


Fig.  4.  Equivalent  electrical  circuit. 
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where  Z)eff  and  cg  are  the  effective  gas  diffusion  coefficient  and 
the  molar  concentration  of  gas  mixture  defined  respectively  as, 


o'ff  =  £*(1  -  s)L5d‘J 


(9) 


Cg  ~ 


(10) 


This  gas  pressure  of  Pg  can  be  obtained  from  the  multiphase 
mixture  model  [4],  as  rather  difficult  function  of  saturation  of  s 
through  capillary  pressure: 


ps 


kd 

V] 


dPCA 

v - as 

ds 


(ID 


where  kY  is  phase  change  rate  [11].  This  expression  means  that 
the  phase  change  is  progressed  within  finite  rate  and  that  it  allows 
supersaturation  state. 

3.2.4.  Equivalent  electric  circuit 

We  considered  the  equivalent  electric  circuit  [10]  shown  in 
Fig.  4.  The  circuit  elements  are  coupled  with  the  gas  species  con¬ 
centrations  in  GDL  and  with  the  water  content  in  PEM,  which  are 
obtained  by  the  conservation  equations  mentioned  above.  Solv¬ 
ing  this  circuit  by  Kirchhoff ’s  law,  the  distribution  of  cell  voltage, 
overpotential  and  current  density  distribution  are  obtained. 

Local  cell  voltage  is  calculated  from  the  following  equation. 


3.2.2.  Supplemental  relations  and  mixture  properties  [4] 

In  GDL  of  porous  media,  the  gas  and  the  liquid  phase  veloc¬ 
ities  denoted  as  iq  and  ug  are  expressed  in  the  following. 


PgUg  —  —  ji  +  Ag/ou 
pm  =  ji  +  ^ipu 


In  GFC,  we  assume  that  the  two-phase  flow  is  well  mixed:  the 
gas  phase  velocity  is  equal  to  liquid  phase  velocity.  The  capillary 
flux  ji  can  be  expressed  as, 


Ji  = 


AiAg  K 


VPC 


The  capillary  pressure  Pc  is  given  by, 


/  £  \  0-5 

Pc  =  pg  —  P|  =  ^  —  J  a  <zos6cJ(s) 


(13) 


(14) 


The  Leverett’s  function  of  J(s)  takes  the  following  form: 


J  1.417(1  -  s)  -  2.12(1  -  s)2  +  1.263(1  -  s)3, 
\  1.417s  —  2.12s2  +  1.263s, 


0C  <  90° 
6>c  >  90° 


(15) 


The  coefficients  of  mixture  transport  properties  appeared  in 
above  equations  are  given  as: 


Ecell  —  Voc  I  Rohm  +  ?7cat  (20) 

Here,  Voc  is  the  open  circuit  voltage  and  ?7cat  is  the  cathode 
activation  overpotential.  Voc  is  calculated  from  following  Nernst 
equation. 

Voc  =  1-229  -  8.456  x  10_4(7,ceu  -  298.15) 

/  p  p0.5  \ 

+  4.308  X  10-5rceii  In  °2  (21) 

V  PH20  J 


Rohm  is  the  Ohmic  resistance,  and  it  mainly  consists  of  proton 
resistance  in  PEM.  Therefore  R0hm  can  be  obtained  as  the  func¬ 
tion  of  local  proton  electric  conductivity,  namely  as  the  function 
of  local  water  content  in  PEM: 

r8m  j 

^ohm  —  /  s  w\  A  ^  (22) 

JO  °rm(Pm)Alact 

In  equivalent  electric  circuit  as  shown  in  Fig.  4,  the  CL  is 
denoted  as  a  parallel  circuit,  which  consists  of  electric  double 
layer  capacitor  and  nonlinear  activation  overpotential  resistance. 
The  electric  charge  accumulated  at  the  capacitor  is  calculated 
from, 


Relative  permeability: 
kx i  =  s3,  kxg  =  (1  —  s)3 
Kinematic  viscosity: 

1 

v  =  - 

(Wh)  +  kvg/vg 

Relative  mobility: 

kY\ 

Ai  =  — v,  Ag  =  1  —  Ai 


(16) 


(17) 


(18) 


3.2.3.  Phase  change  rate 

The  phase  change  rate  mgi  in  Eqs.  (5)  and  (6)  can  be  expressed 
as: 


d_Q 

dt 


=  (i  /BvMact 


(23) 


where  i  is  the  local  current  density.  The  /bv  in  this  equation  is 
the  current  through  the  nonlinear  resistance,  and  it  is  calculated 
from  a  Butler- Volmer  equation. 


z’bv  =  (1  -  s)2ioc 


exp 


(24) 


In  steady  state,  the  local  current  density  corresponds  to  this 
current  density  of  /bv- 

Cathode  activation  overpotential  r]CSLt  is  given  as: 


7?  cat  — 


Gcat 

Ccat 


(25) 


OTgi  =  fcv(pvap  -  Psa?)/^)’ 

'  l  -  S 
s. 


m  = 


if  (Pvap  -  PZP)  >  0 

if(/0Vap  -  Psat)  <  0 


where  Ccat  is  the  capacitance  of  the  electric  double  layer  capac¬ 
itor,  which  changes  with  water  content  at  CL.  Anode  activation 
overpotential  is  neglected  here,  because  it  is  enough  small  not 
to  influence  the  current  distribution. 
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3.3.  Boundary  conditions 


At  the  inlet  of  GFC,  constant  concentration  of  each  species 
and  of  flow  rate  is  specified.  The  inlet  velocity  of  GFC  is  defined 
as  a  function  of  representative  current  of  /rep.  The  inlet  velocity 
for  anode  and  cathode  GFC  are  respectively, 


u  in,  cat 


/rep/(4F) 

§catT^2^gfc 


win,ano 


Wg  F) 

§ano^2^gfc 


(26) 


Where  Ac ha  is  the  channel  cross-section  and  §  is  the  utilization 
ratio  of  reactant  gases  to  the  representative  current. 

At  the  outlet  of  GFC,  it  is  assumed  that  the  velocity,  the 
concentration  and  the  liquid  saturation  fields  are  fully  developed, 
and  their  gradients  are  zero: 


3u 

3  cl 

ds 

—  =0, 

—  =0, 

3  n 

3  n 

3  n 

(27) 


where  n  is  normal  vector  to  the  boundary  surface.  At  the  bound¬ 
ary  between  separator  and  GFC,  non-slip  and  impermeable 
conditions  are  applied  for  the  mixture  of  gas  and  liquid. 


u  =  0, 


(28) 


At  the  boundary  between  GDL  and  PEM,  non- slip  and  imper¬ 
meable  conditions  are  also  applied  for  the  mixture  of  gas  and 
liquid. 


u  =  0, 


(29) 


This  pw|vap,eq  agrees  with  the  value  reported  by  Springer  et 
al.  [14]  within  20%,  meaning  that  Eq.  (31)  roughly  covers  the 
equilibrium  state. 

We  also  well  consider  the  liquid  water  transport  through 
GDL/GFC  interface  with  the  modification  of  the  transport  model 
driven  by  capillary  pressure  [11].  Liquid  water  flux  from  GDL 
to  GFC  is  given  as 

jib  =  -«bpliq  —  (  -  VVs  (32) 

Ml  V  <W 

In  this  equation,  we  added  a  variable  of  n\>  and  considered  flow 
rate  dependence  of  jib  with  this  variable.  With  this  modification, 
the  numerical  simulation  in  the  present  study  seems  to  show 
the  better  agreement  with  measurement  result.  The  derivative  of 
(— dPc/ds)  is  treated  as  a  constant  value  following  the  previous 
study  [11],  though  it  is  given  as  the  function  of  saturation  of  s 
in  general. 

3.5.  Numerical  procedure 

All  equations  shown  in  this  section,  which  are  mass  and 
momentum  conservation  equations,  equivalent  electric  circuit, 
and  supplemental  relations,  are  coupled  and  solved  numeri¬ 
cally.  The  equations  of  continuity  and  momentum  conservation 
equation  are  solved  by  the  SIMPLE  method  [15].  Physical,  elec¬ 
trochemical  and  transport  properties  for  simulation  are  listed  in 
Table  4. 

4.  Experimental  results  and  discussion 


3.4.  Water  transport  model  at  interface 

To  obtain  father  agreement  between  numerical  solution  and 
visualization  result,  we  especially  pay  attention  to  the  water 
transport  at  the  PEM/GDL  and  the  GDL/GFC  interface.  We 
consider  the  non-equilibrium  water  content  in  PEM  near  the 
PEM/GDL  interface,  and  distinguish  liquid  water  transport  from 
vapor  transport.  Including  these  concerns,  the  expression  of 
water  flux  is  given  with  a  modification  of  previous  study  of 
[12]. 

fw  =  (1  —  ‘S’)^vapPgCQ  —  (1  —  ^)^m,vapPm  T"  ^liqPl  —  ^m,liqPm 

(30) 


The  first  term  in  RHS  is  the  water  vapor  flux  through 
PEM/GDL  interface,  which  starts  from  the  vapor  phase  in  GDL. 
The  second  term  is  the  counter  flux  to  the  first  term.  The  third 
term  and  the  forth  term  are  liquid  water  flux  from  GDL  and  from 
PEM,  respectively.  The  vapor  flux  coefficients  of  kVap  and  km ;  vap 
were  given  in  [12].  The  liquid  flux  coefficients  of  k\[q  and  &m,liq 
were  estimated  from  Zawodzinski  et  al.  [13].  The  equilibrium 
water  densities  in  PEM,  in  the  cases  that  the  GDL  is  filled  with 
water  vapor  and  with  liquid  water,  are  easily  derived  from  Eq. 
(30): 


Pm  I  vap,  eq  — 


'Wap  vap 

t — pJ> 

Vm,vap 


Pmlliq,eq  — 


_^liq_pliq 

^m,liq 


(31) 


In  this  section,  we  show  cell  voltage,  and  the  visualized  image 
of  flooding  in  cathode  GFC,  and  their  time  evolution.  Based  on 
these  experimental  results,  we  discuss  the  influence  of  operation 
condition  on  flooding  and  on  cell  voltage. 

4.1.  Influence  of  flow  rate  on  flooding 

Fig.  5  is  time  evolutions  of  cell  voltage,  under  the  same  cur¬ 
rent  density  of  0.6  A  cm-2,  with  different  flow  rate  of  58,  116 
and  232  seem,  corresponding  to  the  experimental  condition  of 
No.  3-1,  3-2  and  3-3  as  shown  in  Table  2.  Fig.  6  is  the  visu¬ 
alized  images,  which  are  taken  during  the  cell  operation  under 
these  three  conditions.  These  figures,  as  a  whole,  indicate  that 
the  vapor  pressure  in  the  cell  reached  saturation  pressure,  and 
flooding  area  gradually  spread  in  GFC,  resulting  in  the  reduction 
of  cell  voltage. 

At  the  moment  when  the  current  was  loaded  from  0  to 
0.6  A  cm-2,  the  cell  voltage  instantly  dropped  as  shown  in  Fig.  5, 
and  this  instant  voltage  drop  was  enhanced  by  smaller  flow 
rate.  This  flow  rate  dependence  on  the  cell  voltage  might  be 
caused  by  the  shortage  of  oxygen:  smaller  flow  rate  leads  to  the 
shortage  of  oxygen  at  cathode  CL,  resulting  in  larger  activation 
overpotential. 

After  the  instant  voltage  drop,  each  trend  of  cell  voltage 
showed  relatively- slow  decrease,  and  smaller  flow  rate  invited 
higher  decreasing  rate.  Smaller  flow  rate  seems  to  deteriorate  the 
drainage  of  liquid  water  in  GDL  and  GFC,  leading  to  the  high 
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Table  4 

Physical,  electrochemical  and  transport  properties 
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Description 

Value 

Reference 

Anode/cathode  pressure,  P  (atm) 

1. 0/1.0 

Contact  angle  of  GDL,  6C  (°) 

110 

Electro-osmotic  coefficient,  n d 

1.0 

Faraday  constant,  F  (C  mol-1) 

96487 

Liquid  water  density,  p\  (kgm-3) 

995.7 

Permeability  of  GDL,  K  (m2) 

1.0  x  10-12 

Surface  tension,  o  (N  m-1) 

0.0627 

Universal  gas  constant,  R  (J  (molK)-1) 

8.314 

Viscosity  of  liquid  water,  fi\  (Pas) 

0.8  x  10-3 

Water  diffusivity  in  PEM,  Dm  (m2  s-1) 

1.5  x  10-6 

Cathode  exchange  current  density,  z'oc  (A  cm-2) 

0.0017 

Measured 

Liquid  water  transport  coefficient  between  GDL/GFC,  n\> 

0.08  (Air  232  seem) 

0.04  (Air  116  seem) 

0.02  (Air  58  seem) 

Adjustable 

—dPc/ds,  (Pa) 

22.95 

Liquid  water  permeability,  K\  (m2) 

1.1  X  lO"13 

[11] 

Phase  change  rate  constant,  kv  (1  s-1) 

100 

Vapor  transport  coefficient  from  PEM  to  GDL,  km?va p  (ms-1) 

1.46  x  10“7 

[12] 

Vapor  transport  coefficient  from  GDL  to  PEM,  kvap  (ms-1) 

1.58  x  10“3 

Liquid  transport  coefficient  from  PEM  to  GDL,  km, nq  (ms-1) 

8.24  x  10“6 

Estimated  from  [13] 

Liquid  transport  coefficient  from  GDL  to  PEM,  k\iq  (ms-1) 

6.01  x  10-6 

Effective  porosity  of  GDL,  £*(=£/ r) 

0.21  (cathode) 

0.09  (anode) 

[16] 

decreasing  rate  of  cell  voltage.  In  addition  an  immediate  voltage 
change,  which  is  described  as  a  vertical  line  in  Fig.  5,  was  also 
observed  especially  under  smaller  flow  rate  conditions  (No.  3- 
3,  3-2).  Thus  the  trends  of  cell  voltage  in  our  experiment  were 
categorized  into  two  types:  the  relatively- slow  decrease  and  the 
immediate  change. 


Fig.  5.  Time  variations  of  the  cell  voltage  under  the  current  condition  of 
0.6  A  cm-2,  with  the  different  air  flow  rate:  232  seem  (No.  3-3):  1 16  seem  (No. 
3-2):  58  seem  (No.  3-1). 


Fig.  6(a-d)  are  the  successive  image  of  water  behavior  in 
cathode  GFC  under  0.6  A  cm-2  and  232  seem  condition  (No.  3- 
3).  After  5  min  elapsed  from  applying  load  current,  small  water 
droplets  began  to  appear  near  the  outlet  of  GFC  as  shown  in 
Fig.  6a,  and  the  droplets  almost  popped  up  from  the  edge  of 
GFC.  Liquid  water  is  thought  to  prefer  to  accumulate  under  the 
rib,  where  less  convection  is  expected.  It  is  noted  that,  despite 
experimental  conditions,  the  point  where  water  droplets  pop  up 
was  exactly  same.  This  means  that  the  dominant  path  to  flow 
liquid  water  exists  in  GDL,  and  that  liquid  water  selectively 
migrates  there. 

After  10  min  and  30  min  elapsed,  the  water  droplets  grew, 
and  combined  near  the  outlet  and  at  the  midstream  of  GFC. 
In  addition,  the  surface  of  glass  was  fogged  partially.  This 
fogged  area  might  correspond  to  supersaturated  area.  After  tak¬ 
ing  Fig.  6c,  accumulated  water  flushed  out  of  GFC,  and  the 
cell  voltage  immediately  recovered  by  about  40  mV.  Fig.  6d  was 
taken  just  after  this  event,  and  it  was  found  that  the  water  droplets 
were  cleaned  off  by  the  flushing.  The  above  step-by-step  events, 
namely  the  droplets  popping-up,  their  growth,  their  combination 
and  flushing,  were  observed  repetitively,  and  it  is  thought  to  be 
one  of  the  characteristic  phenomena  under  a  flooding  condition. 

Fig.  6e  and  f  are  the  images  under  the  0.6  A  cm-2  and 
116  seem  condition  (No.  3-2).  Fig.  6e  shows  that  water  droplet 
growth  under  this  condition  was  faster  than  that  under  232  seem 
condition  (No.  3-1).  Comparing  between  Fig.  6c  and  f,  it  can 
be  seen  that  the  more  water  in  case  of  slower  flow  rate  was 
accumulated  in  GFC  even  in  shorter  period.  Just  after  Fig.  6f 
was  taken,  the  accumulated  water  flushed.  Different  from  the 
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Fig.  6.  Images  of  cathode  GFC  at  the  specified  time,  under  the  current  condition  of  0.6  A  cm  2,  with  the  different  air  flow  rate:  232  seem  (No.  3-3):  1 16  seem  (No. 
3-2):  58  seem  (No.  3-1). 


232  seem  condition,  cell  voltage  immediately  dropped,  and  elec¬ 
tricity  generation  stopped  after  the  flush.  When  this  flushing 
happened,  it  was  observed  that  something  like  water  film  cov¬ 
ered  the  large  area  of  GDL  surface.  It  can  be  thought  that  this  film 
triggered  by  the  flushing  fully  clogged  the  pathway  of  oxygen 
gas  supply  to  CL  and  invited  the  immediate  voltage  drop  result¬ 
ing  in  the  electricity  generated  stop.  The  liquid  film,  mentioned 
just  above,  formed  and  disappeared  within  1  s,  and  it  emerged 
at  irregular  intervals.  On  the  other  hand  the  minimum  intervals 
between  recordings  image  is  15  s  with  the  microscope  used.  This 
restriction  invited  the  difficulty  to  record  the  image  of  liquid  film. 
Regrettably  we  could  not  succeed  to  record  the  image  of  liquid 
film,  although  we  caught  it  visually. 

Fig.  6g  is  the  image  under  the  0.6  A  cm-2  and  58  seem  con¬ 
dition  (No.  3-1).  Under  this  smallest  flow  rate  condition,  the 
voltage  decreasing  rate  was  highest.  After  2  min  of  operation, 
cell  voltage  decreased  below  0.15  V  and  we  stopped  the  opera¬ 
tion.  During  this  2  min,  no  droplet  was  observed  in  GFC.  This 
voltage  drop  did  not  seem  to  be  caused  by  only  a  shortage  of 
oxygen,  because  gas  re-distribution  after  applying  load  current  is 
estimated  as  to  need  only  about  0.01  s.  This  representative  time 
supports  the  explanation  that  the  flooding  in  GDL  and  CL,  which 
takes  several  minutes  to  develop,  caused  the  voltage  decreasing. 

Summing  up,  it  found  that  relatively- slow  voltage  decrease 
and  the  immediately  voltage  change  are  respectively  caused  by 
the  flooding  in  porous  media  (GDL  and  CL)  and  by  the  flushing 
in  GFC. 

4.2.  Influence  of  current  on  flooding 

In  this  sub-section,  the  influence  of  load  current  on  flooding  is 
discussed  by  the  experimental  results  of  Figs.  7  and  8,  under  the 


same  flow  rate  condition  of  116  seem,  with  the  current  density 
of  0.15,  0.3  and  0.6  A  cm-2,  corresponding  to  the  condition  of 
No.  1-2,  2-2  and  3-2  shown  in  Table  2. 

As  a  whole,  it  is  found  that  a  larger  load  current  invited  the 
wider  flooding  area  and  the  higher  decreasing  rate  of  cell  volt¬ 
age.  In  the  case  of  minimum  current  condition  of  0.15  A  cm-2, 
the  cell  voltage  almost  kept  constant,  and  the  flooding  in  GFC 
was  not  developed.  In  the  cases  of  larger  current  condition  of  0.3 
and  0.6  A  cm-2,  the  cell  voltage  kept  decreasing,  and  the  flood¬ 
ing  was  developed  gradually.  After  the  14  and  36  min  elapsed  in 


Elapsed  time  /  min 


Fig.  7.  Time  variations  of  the  cell  voltage  under  the  air  flow  rate  condition  of 
116  seem,  with  the  different  current:  0.6  A  cm-2  (No.  3-2):  0.3  A  cm-2  (No. 
2-2):  0.15  Acm“2  (No.  1-2). 
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Fig.  8.  Images  of  cathode  GFC  at  the  specified  time,  under  the  air  flow  rate  condition  of  116  seem,  with  the  different  current:  0.3  A  cm  2  (No.  2-2):  0.15  A  cm  2 
(No.  1-2). 
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Fig.  9.  Time  variations  of  the  cell  voltage  under  the  utilization  ratio  condition 
of  0.05,  with  the  different  current  and  air  flow  rate. 


and  flow  rate,  whose  experimental  numbers  are  No.  1-1, 2-2,  3-3 
as  shown  in  Table  2.  Contrary  to  our  expectation,  Figs.  9  and  10 
indicate  that  the  cell  voltage  under  the  highest  current  (largest 
flow  rate)  condition  is  stable  rather  than  that  under  the  lower 
current  (smaller  flow  rate)  condition.  This  could  be  caused  by 
the  effect  that  large  flow  rate  enhance  the  drainage  of  water 
droplet  in  GFC.  This  flow  rate  effect  is  confirmed  by  the  image 
of  flooding  in  Fig.  10.  Comparing  the  image  between  different 
flow  conditions,  the  flooding  area  in  case  of  the  highest  current 
(largest  flow  rate)  condition  is  small,  though  the  utilization  ratio 
is  same.  Thus  flow  rate  is  one  of  key  parameter  to  determine  the 
degree  of  flooding. 

The  immediate  voltage  drops,  shown  in  the  case  of  No.  1- 
1  and  2-2  in  Fig.  9,  happened  again  after  the  flushing,  being 
same  fashion  as  shown  in  previous  sections.  After  load  current 
applied,  water  droplets  popped  up  in  GFC,  and  they  accumulated 
there,  and  were  flushed,  resulting  in  the  large  voltage  drop.  On 
the  other  hand  the  cell  voltage  was  rather  stable  in  case  of  largest 
flow  rate  of  No.  3-3.  In  this  case  water  seems  to  be  well  balanced 
with  repetitive  accumulation  and  drainage. 


cases  of  these  higher  0.3  and  0.6  A  cm-2  condition  respectively, 
the  cell  voltage  dropped  immediately  and  the  electricity  genera¬ 
tion  stopped  with  the  flushing  as  also  mentioned  in  Section  4.1, 
though  the  flooding  did  not  fully  covered  the  GDL  surface  as 
shown  in  Fig.  8.  Thus  it  is  understood  that  the  flushing  is  a  key 
event  to  cause  the  immediate  generation  stop. 

4.3.  Comparison  of  flooding  between  same  utilization  ratio 
conditions 

In  this  sub-section,  the  flooding  are  compared  between  same 
utilization  ratio  conditions  of  0.05  with  changing  load  current 


5.  Numerical  simulation  results  and  its  comparison  with 
experimental  results 

5.1.  Time  evolution  of  cell  voltage 

Fig.  1 1  shows  the  cell  voltage  change  for  nine  cases  of  oper¬ 
ation  condition  as  listed  in  Table  2.  Fig.  11a  and  b  are  the 
experimental  and  the  simulation  result,  respectively.  Different 
current  density  condition  is  indicated  as  dotted  circle  in  the  fig¬ 
ures.  For  each  current  density  condition,  three  voltage  trends 
corresponding  to  three  flow  rate  conditions  are  plotted.  As  a 
whole  of  the  experimental  result,  the  cell  voltage  was  high  under 


0.6 
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Fig.  10.  Images  of  cathode  GFC  at  the  specified  time,  under  the  utilization  ratio  condition  of  0.05,  with  the  different  current  and  air  flow  rate. 
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Fig.  11.  Comparison  of  the  cell  voltage  change  between  experiment  and  simulation,  (a)  Experimental  results,  (b)  Simulation  results. 


lower  current  density  and  higher  flow  rate,  although  the  flow 
rate  dependence  was  rather  obscure  in  cases  of  small  current 
condition.  Our  simulation  well  predicts  this  overall  trend  of  cell 
voltage  change. 

In  case  of  0.15  A  cm-2  of  the  experimental  result,  cell  volt¬ 
age  kept  high,  and  flow  rate  dependence  of  it  was  small.  The 
simulation  result  agreed  with  the  experimental  result. 

In  case  0.3  A  cm-2,  experimental  results  showed  that  the  cell 
voltage  slightly  reduced  as  time  elapsed,  and  that  its  reduction 
rate  increased  as  the  air  flow  rate  decreased.  The  simulation 
result  seems  to  replicate  these  features.  However,  the  simulation 
result  could  not  replicate  the  immediate  voltage  drop,  which 
emerged  in  case  of  58  and  116  seem  flow  rate.  As  mentioned  in 
the  previous  section,  this  immediate  voltage  drop  was  trigged  by 
the  flushing  event,  which  successively  happens  after  the  flooding 
develops  in  GFC.  On  the  other  hand,  the  simulation  is  not  able 
to  simulate  the  flushing  event,  because  gases  and  liquid  water 
in  GFC  well  mixed  and  flow  together  without  flooding  in  GFC. 
Thus  simulation  is  not  able  to  replicate  the  immediate  voltage 
drop. 

In  case  of  0.6  A  cm-2,  experimental  result  showed  that  the 
dependence  of  the  flow  rate  on  cell  voltage  was  large,  and  that 
the  cell  voltage  decreased  rather  quickly  within  15  min  in  the 
lower  flow  rate  cases.  The  simulation  result  replicated  these  fea¬ 
tures  qualitatively.  Fig.  12  obtained  by  the  simulation  shows 
the  liquid  water  saturation  and  oxygen  gas  density  at  cathode 
CL.  In  the  case  of  highest  flow  rate  of  232  seem,  water  bal¬ 
ance  between  accumulation  and  drainage  seems  to  be  acquired, 
resulting  that  the  liquid  saturation  and  the  oxygen  gas  density 
reached  constant  value.  On  the  other  hand,  the  water  in  the  case 
of  smaller  flow  rate  did  not  balanced,  and  then  the  liquid  satu¬ 
ration  and  the  oxygen  gas  density  at  cathode  CL  continuously 
increased  and  decreased,  respectively.  Namely  the  flooding  in 
cathode  GDL,  especially  near  cathode  CL,  developed,  and  liq¬ 
uid  water  droplets  hindered  the  oxygen  transport  in  GDL,  and 
then  the  cathode  activation  overpotential  increased,  resulting  in 
the  quick  decrease  of  cell  voltage  shown  in  the  experimental 
result. 


Putting  together  the  above  consideration,  it  is  understood  that 
the  immediate  voltage  change  is  triggered  by  flushing,  and  that 
it  is  not  able  to  replicate  with  present  our  simulation.  On  the 
other  hand,  it  is  found  that  the  relatively- slow  voltage  drop  is 
caused  by  the  flooding  growth  in  GDL  and  CL,  and  that  it  can 
well  replicate  with  our  simulation. 

5.2.  Comparison  of  GFC  image  between  simulation  and 
experiment 

Fig.  13  is  the  comparison  of  GFC  images  between  the  exper¬ 
iment  and  the  simulation.  Each  image  was  taken  at  the  time 
specified  in  Fig.  11.  The  simulation  image  expresses  vapor 
density  distribution  with  the  highlighting  where  vapor  is  super¬ 
saturated  or  unsaturated.  In  the  experiment  image,  the  liquid 
droplet  and  fog  are  spread  from  the  outlet  to  the  upstream  in 
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Fig.  12.  Time  variation  of  the  liquid  water  saturation  and  the  oxygen  gas  density 
at  the  cathode  CL  (the  interface  between  PEM  and  cathode  GDL)  in  case  of 
0.6  A  cm-2. 
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Fig.  13.  Comparison  between  the  cathode  image  and  the  vapor  density  distribution  obtained  by  visualization  experiment  and  numerical  simulation,  respectively. 


GFC,  and  the  region  becomes  larger  as  the  air  flow  rate  decreases. 
The  simulation  result  seems  to  replicate  these  features,  espe¬ 
cially  on  the  beginning  point  of  the  supersaturation  region.  Thus, 
it  was  confirmed  that  the  observed  fogging  area  corresponded 
to  the  supersaturated  area  in  GFC. 

6.  Conclusion 

We  developed  a  two-phase  flow  and  2D  simulator,  where 
water  transport  at  the  interface  of  PEM/GDL  and  GDL/GFC 
was  carefully  treated.  Numerical  simulation  was  executed  and 
compared  with  visualization  experiment,  reading  to  the  follow¬ 
ing  results.  The  Simulation  could  replicate  the  dependence  of 
load  current  and  flow  rate  on  cell  voltage.  The  simulation  could 
also  replicate  the  relatively-slow  decrease  of  cell  voltage,  which 
was  caused  by  the  flooding  in  GDL  and  CL.  Moreover,  the  sim¬ 
ulation  could  replicate  the  supersaturation  region  in  GFC  and 
the  flow-rate  dependence  of  it. 

Visualization  experiment  caught  the  immediate  cell  voltage 
change  with  liquid  water  film,  which  momentarily  formed  and 
covered  GDL  surface,  synchronizing  with  the  flushing-out  of  liq¬ 
uid  water  accumulated  in  GFC.  This  immediate  voltage  change 
with  the  dynamic  water  behavior  was  thought  to  be  one  of 
representative  phenomena  under  flooding  condition.  Our  sim¬ 
ulation  could  not  replicate  this  phenomenon,  because  of  the  use 
of  well-mixed  modeling  in  GFC. 
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